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We study the dynamics involved in a sparse random network model. We extend the standard 
mean-field approximation for the dynamics of a random network by employing the path-integral 
approach. The result indicates that the distribution of the variable is essentially identical to that 
obtained from globally coupled oscillators with random Gaussian interaction. We present the results 
of a numerical simulation of the Kuramoto transition in a random network, which is found to be 
consistent with this analysis. 



I. INTRODUCTION 

a . 

Many systems in nature, such as food webs, metabolic systems, coauthorship of papers, the worldwide web, and so 
C/3 ' on, can be represented as complex networks 0,0,0,0. Investigations of real networks have shown that these networks 
have topologies different from random networks. In particular, we have recognized that many networks have scale-free 
degree distribution, P{k) cx fc~ 7 , where k is the degree of nodes. 

The dynamics involved in complex networks has become an important aspect of the complex network studies in 
recent times. This problem includes, for example, the spreading of virus in the internet, synchronization of neurons in 
a brain, change of populations in a food web. Recently, Pastor-Satorras and Vespignani obtained unexpected results 
in this regard ja] • They studied the spread of viruses in complex networks and found that no threshold of infection 
rate exists for the susceptible- infected-susceptible model in the random scale-free network with 7 < 3, if the size of 
the network N is infinite. Though real networks such as the internet are finite-size network, this result implies that a 
virus with a small infection rate can spread over the whole network. We had previously presented another remarkable 
example of the unusual dynamics involved in complex networks0,0- We studied the Kuramoto synchronization in 
a random network of oscillators and found that the critical coupling for synchronization becomes zero in scale-free 
network with 7 < 3. 

In these studies, the mean-field approximation plays an essential role. For the mean-field approximation, we consider 
' a model in which a node i couples to another node j with a strength proportional to "mean coupling probability" 
, hkj/ktot., where ki and k to t. are the degree of node i and total number of edges, respectively. The dynamics in 
complex network is much simplified by this approximation, and we can obtain analytical results. However, this model 
differs from the original network model, in which each node couples to a finite number of nodes. It is remarkable that 
the mean-field approximation is in good agreement with the numerical simulation result of a random network model. 

One of the objectives of this paper is to provide a sound explanation for the mean-field approximation. It is 
unclear why the mean-field approximation performs well in the random network model. The validity of the mean-field 
approximation, particularly with regard to the Kuramoto transition, is debatable. Moreno, Pacheko and Vazquez- 
Prada carried out numerical simulations on the Kuramoto transition in the Barabashi- Albert network. 0. They 
concluded that the critical coupling K c is not even if N — > 00. Their conclusion seems to contradict the result of 
mean-field theory K c = 0, though this discrepancy is possibly due to the difference of the order parameter used in 
O ■ these papers. Restrepo, Hunt, and Ott suggested that the the argument based on the largest eigenvalue of the network 
matrix is superior to that based on the mean-field theory They demonstrated that the mean-field approximation 
is valid for Erdos-Renyi networks and random scale-free networks with 7 = 3, while this approximation does not hold 
in the case of scale-free networks with 7 = 2. However, they did not provide any explanation as to why the mean-field 
theory works well in some random network models. An appropriate explanation to this question is a matter of great 
interest and significance. 

The second objective of this study is to extend the mean-field theory. Although the mean-field theory displays good 
qualitative coincidence with numerical simulation, it is impossible to examine the fluctuation of the variables by the 
mean-field approximation. Moereover, as noticed, we cannot apply the mean-field approximation in some network 
models. Therefore it is meaningful to make an approximation that covers a wider range of complex networks. In this 
paper, we demonstrate that the distribution of variables in the sparse random network model can be approximated by 
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that obtained from a globally coupled network, in which the distribution of the interaction between the nodes is given 
by a Gaussian random number. This result indicates that the dynamics in random network can be approximated 
more precisely by appropriate methods such as dynamical mean- field theory [13 . 

In order to realize the above-mentioned objectives, we utilize the path-integral approach. The path integral, 
which was originally developed for application in quantum mechani cs 1 12). h as also been applied to random impurity 
problems^, Q], random spin glasses^, 0, 0], neural networks|l8l llflj. and oscillator systems j^- One of the 
advantages of this approach is that the average over an ensemble of networks can be calculated easily. Limitations of 
the path-integral include an infinite number of integrals and obtaining a precise average over the ensemble, which is 
not usually possible. However, this method enables approximation of the distribution of the variables in a systematic 
manner. In particular, the mean-field approximation can be derived as the lowest order approximation of the path 
integral. The methods used by us are similar to that used by Theumann for the Hopfield network model [22|. 

The outline of this paper is as follows. In Sec. [HI we present the general description of the dynamics of a 
network model based on the path-integral approach. We derive a formula that is general and can be applied to any 
network model in this section. In Sec. IIHI we present two approximations of the path-integral formula, mean-field 
approximation and perturbation. We also prove that the dynamics of a random network is essentially identical to 
that of a random Gaussian network. In Sec. IIVI we apply the analysis to the Kuramoto transition in a random sparse 
network. We present the results of numerical simulation, which are consistent with that obtained from the analysis. 
To conclude, we discuss our obtained results. 



II. PATH-INTEGRAL APPROACH TO THE DYNAMICS OF A NETWORK MODEL 

In this section, we introduce the formalism to study the dynamics of a network model using the path-integral 
approach. We consider the following differential equations for the network model: 



A" 



Xi = fjjxj) + y] ai,jg(xi, Xj) + &(t), 

3=1 



(1) 



where £i(t) is a random force that satisfies = 0), (£i(t)£j(t )) = Si,jS(t — t )a 2 . We assume Xi — ai^o at t = 0. In 

order to discuss the dynamics of this system, it is useful to introduce Matrin-Siggia-Rose(MSR) generating functional 
Z, which is defined as[HIT^| 



/ j \ NN t I /. N N t \ 

Z[{k tk }, {k,k}] = f - ] ( / II II dx i,kdx iik e~ s exp(l i: kXi ik + h,kXi,k)J 

" 1 W i=l k=0 I 



where the action S is given by 
a 2 At . 



k + ix itk I x iik - a^fc-i - At(fi(x it k-i) + ^ ai,jg{xi,k-i,%j,k-i)) 



and (• • •) represents the average over the ensemble of networks. J is the functional Jacobian term, 



(2) 



(3) 



J = exp I 



d{fi{xi,k) + a i::j g(x i:k ,x jt k)) 



Ox 



i , k 



(4) 



Though this term is necessary for the renormalization Z(0) — 1, it is a little cumbersome to treat it in a practical 
calculation, such as the mean-field approximation or a perturbation. Here we note that, as De Dominicis showed [T^. 
the only effect of this Jacobian term is to subtract the nonretarded correlation function (xi^Xj^k+k'}, where k' > 1. 
In the following discussion, we omit this Jacobian term, remembering that we only consider the retarded correlation 
function. Maintaining AtN t constant at the limit At — > 0, we obtain the MSR generating functional. 
We consider the network described by 



1 with probability pi 
with probability 1 — p 



(5) 



3 



We note that py can be a function of variables such as i or j. For example, in the one-dimensional chain model, py 
is 1 if \i — j\ = 1, else it is 0. The average over all networks can be expressed as 



exp 



i , k 



}=n 


Pij exp | 


1 1,3 





i j exp iAtoj,fcg(zj,fc_i, gj.fc-i) j + 1 -py 



(6) 



and we obtain 



where 



exp( 



S )H 



h3 



Pij exp |y^iAta;,- l jfcff(a;t ) fc-i,a;j,A_i)| + 1 



x a 2 At 

i.k 



Ik + i%i,k{xi,k - Xi,k-i - At/i(aFj,fc_i)}. 



(7) 



(8) 



The above-mentioned expression is a general one and can be applied to the dynamics of any network model. However, 
it is often impossible to calculate the precise value of (e~ s ), particularly in the case of nonlinear dynamics. We need 
an approximation to obtain the value of (e~ s ). In the next section, we approximate Eq. © by assuming p^. <C 1 
and Pi_jPk,i -C Pij for any i,j, fc, and I. 



III. APPROXIMATION OF THE MSR GENERATING FUNCTIONAL IN A SPARSE RANDOM 

NETWORK MODEL 



In this section, we develop an approximation for the MSR generating functional Z in a sparse random network. 
For this, we assume py <C 1 and PyPfe,; <C py for any k, and I. In the case of the Erdos-Renyi model, py is 
independent of i and j; py = (k)/N. Therefore this assumption is valid for a sparse Erdos-Renyi model, because 
Pi,jPk,i = (k) 2 /N 2 <C Py. In the case of random network with distribution P(k), we construct the network as follows. 
First, we define the "degee" of node i as k[, whose distribution concides with P{k). Second, we connect the nodes 
i and j with probability py = k^k'j / ^2 i k' i . Using this procedure, we obtain the random network whose degree 
distribution is approximately given by P(k). In this case, if the maximum degree of a node k max is much smaller than 
N, k max -C N, the assumption is satisfied. On the other hand, this assumption is not satisfied in the Watts-Strogatz 
model, because p^-j+i ~ 1. 

Since py <C 1, the approximate value of the logarithm of the right-hand side of Eq. © is expressed as follows: 



In ] [ py exp < iAtXj,kg(xi,k-i,Xj,k-i) f + 1 - Pi,j ~ y ~Pi,j + Pij ex P \ iAtXi,kg(xi,k-uXj,k-i) ? 

L K k J J / i,j I k ) 

oo 1 / y 

= y Pi,jy^ 7| I y iAtXj,kg(xi,k-i, Xj,k-i) J • (9) 



Therefore, we obtain 



(e~ s ) ~ exp 



x exp 



E 



i,fc - ix t ,k{xi,k - a?i,fc— i - At/i(a;i ) fc_i)} 



(10) 



To calculate Z from this equation, we need an infinite number of integrals, and we cannot carry out this integration 
practically. However, this formula gives us much information about the averaged dynamics of networks. In the 
following subsections, we consider two simple approximation schemes, the mean-field approximation and perturbation. 
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A. Mean-field approximation and beyond 



To begin with, we consider an approximation that ignores the I > 2 part of Ea. (|10|) and obtain 

a 2 At 



Vexp £ 



{xuk - x^k-i - At(fi(x itk -x) + ^pijgjx^k-i, Xj,k-i))} 



(11) 



This result demonstrates that the MSR generating functional for Eq. JTJ can be approximated as that for the system 
described by 



JV 



(12) 



This equation implies that the mean- field approximation neglects the contribution of the term I > 2 in Ea. (|10|l . 
The mean-field approximation method is based on two assumptions. First, the higher-order term in pij in Eq.@ is 
neglected, and, second, the higher-order term in Eq. (|1U|I is neglected. The former assumption is valid if pij <C 1 for 
all values of i and j. However, neglecting the higher order term is not always valid. In order to examine this argument, 
we study the effect of the term 1 — 2. From the Stratnovich-Hubbard transformation, we obtain the following: 



exp 



Pi.j ~ ( £ iAtx itk g{x lt k-li %j,k-l) 



2npi 



drij exp 



2p t . 



+ i r» 



£ Atx lyk g(x lyk -i- 



Xj,k-l) 



(13) 

By comparing Eqs. (|12ll . and i|13fl . we observe that the MSR generating functional is identical to that of the 

system described by 



JV 



Xi = fi(xi) + ^2(plj + n,j)g{xi, Xj) + £i(t), 



(14) 



where r^j is a random number and its distribution is given by a Gaussian distribution, with a mean value of and a 
dispersion (r 2 j) —Pij- 

Sequential application of the Stratonovich- Hubbard transformation yields the contribution from the term I = 2 n . 
For example, if we consider the term / = 4, then because 



exp 



Pi,j^ (y^JAtx iik g{x iik -i, 



2<7,fc-l, 



2-Kpij 



*i,exp 



2pi,j 



y^jAtxi^gjx^k-i' 



Xj.k—l, 



^ ^ (15) 

the Gaussian fluctuation with dispersion \/pij/3 is added to pi^ in Ea. H13|) . Therefore, (P) is given by the solution 
of 



N 



Xi = fi(xi) + 2j(pi,j + r itj )g(xi, xj) + &(t), 



(16) 



where the distribution of r^j is given by a Gaussian with dispersion a r = J p. 



and the distribution of r' 



is also Gaussian with dispersion a' r — \/pij/3. By sequential application of this transformation, we can obtain the 
effect of the term I = 2™. However, as these correction terms are small for large the Gaussian random network is a 
good approximation of the random sparse network. 

The estimation that utilizes the Stratonovich-Hubbard transformation is very effective in mapping the dynamics of 
the sparse network ensemble onto the dynamics of globally coupled networks. This method is very useful, especially 
in the case where the dynamics of Gaussian random networks is well known. However, the effectiveness of this 
transformation is limited, because it can only consider the terms I — 2™. This method does not elaborate on the 
effect of the terms I = 3,5,6, •• •. In the next subsection we realize that the correction in Z resulting from the term 
I = 2m + 1 is of the order pf j in the random network. 
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B. Perturbation 



As shown in the preceding section, though the mean-field approximation and the Stratonovich-Hubbard transfor- 
mation are very effective methods, they only consider a limited number of terms. To examine the effect of other terns, 
we use the perturbation technique for a network model in this section. The perturbation gives a formal estimate of 
the value of (e~ s ). Furthermore, we note that this method is highly effective because it allows us to estimate the 
accuracy of the approximation by an order of pij. However, this method has two drawbacks. First, it is typically 
impossible to obtain the value of (e~ s ) in nonlinear physics. It is often difficult to obtain the value of (e~ s ) even if 
there is no interaction, and the perturbation can therefore only be applied to limited systems. Second, bifurcation 
or phase transition cannot be obtained without including the infinite order of perturbation in pij. In general, the 
MSR generating functional Z becomes singular at the bifurcation point. However, the finite order perturbation yields 
Z, which is a nonsingular function oipij. Therefore, it is impossible to examine a phase transition by perturbation. 
However, perturbation often yields important information. In this section, we demonstrate that the correction from 
the odd I term is of the order pf j. 

We begin with Eq. H1(J|I . On expanding expQ^pj j " ')> we obtain 

a 2 At 



exp 



E 



oo 1 

y — 



ro=0 



x i,k - iXi,k{Xi,k ~ X ilk -1 ~ Atfifak-l)} 



00 1 ( V 

E Pi J E J\ ( Z^^^hkOi^k-i^j^-i) j 

i,j 1=1 V k ) 



(17) 



We first consider the correction from the term I = 3. Since the effect of the term I = 2 can be expressed as the 
fluctuation in pij , we consider 



exp 



E 

i.k 



fc-i - &t(fi(xi,k-i) + y]p u g(x it k-i , a?j,fc-i))] 



E ~§f ( E i^txi^gixi^- 



l,Xj,k-l) 



(18) 



To account for the effect of the term I — 3, we treat this term using perturbation. We expand the term I — 3 as 



exp 



ly^JAtx itk g(xi t k-i,Xj,k) j = 1 + ^2 -jjf ( ^2iAtxi ik g{x it k- 1 ,Xj t k-i 

i,j \ k / i,j \ k 

+ \ \ E^f (E*^***'^ 35 *'*- 1 ' 3 ^'*- 1 

[ i,j ' V k 




Here we assume Pijp k ,i "C Pij for any value of k, and / again. 
Based on this assumption, (e _s ) can be approximated as 



where 



■At 



1 + E^f (^j^ tx i,k9(xi,k-i,Xj,k-i)^j |e Sl , 



[xi,k - x ijk -i - At(fi(x l:k -i) + y^p ii7 -g(xj,fc_i,a:j-,fe_i))] 



i , k 



To calculate the contributions to Z from these terms, it is convenient to define 



P({ei,k}) 



TJ dx^ k exp ^2 



a 2 At 



ixi,k[xi,k ~ Xi,k-i - M(fi(x itk ^i) 



(19) 



(20) 



(21) 
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(22) 



Since 



„- j. 



-Si 



c) :i 



9ci li k 1 dei l7 k 2 d€i 1 ,k: 



(23) 



the MSR generating functional Z can be calculated from Eq. H2()J) if the differential of P is known. We first consider 
the differential of P in the case where k\, k%, and k^ are distinct. On integrating Eq. I|22|l . we obtain 



P({ei, k })=Cexp 



E 2cr 2At i Xhk ~ Xt ' k ~ i ~ A< I ^(^.fe- 1 ) + EK^ 



£i,fe-l>%',fc-l) — Ci.fc 



(24) 



where C = (f /o^At)^'/ 2 . We note that the integration of P{{e %tk }) over x l)k is 0(f), while P({e iik }) is 0((Ai) 3 / 2 ). 
On differentiating Eq. (|24|l with respect to e^^, we obtain 



^ lx llM - x ll:kl ^! - At \ f n (x nM ^i) +^2p' liJ g(x llM - 1 ,x j ^ 1 ^i)\ \e s \ (25) 



e=0 (7 



where 



E 2^Ai I 



-i - At /j(g,,fc-i) + y^p- i1 -g(xi, fc _i,Xj ! fc_i) 



(26) 



On differentiating Eq. (|25|l with respect to ei lt k 2 and fe 3 we obtain d 3 P / dei ly k 1 dei ly k 2 dei 1 .k 3 . However, we tem- 
porarily consider the term dP/dei^k, because further differentiations with respect to ei 2l k 2 an d £i 3 ,k 3 do not modify 
the following discussion. 

In the limit At — * 0, exp(— S[ ) approaches to the 5 function 5(xi t k—Xi t k—i — At[fi(xi t k-x)+'^2j Pi j9(xi t k-ii %j,k— l)])- 
Therefore, in the limit At — > 0, Eq. Q25|l always attains the value 0. However, Ai7V t is maintained constant when the 
limit At —* is taken. Therefore, the sum of the integrals of dP/dei^k over x^k may have a finite value at the limit 
At — > 0, if Eq.J35J) has a magnitude O(Ai). In order to obtain a more accurate estimate of Eq. (|25fl . we express S[ as 



S 'i ~ E 2(T 2 At i 



+ yfi{xi,k) + y^,Pi,i9( x i,k, x j,k)^j (xi,k+i ~ x i,k) + O(At) + (the terms that do not include x^). (27) 

We consider the integral J dxi^h(xi y k)dP '/ de^k, where h(xi ik ) is an arbitrary nonsingular function. Since the value of 
h^Xi^kjlx^k-Xi^-i-Aty^Xi^-^+p'ijgixi^-iiX^k-i)}] is at x itk = x it k-i + At(fi(x it k-i) +p' iij g{x i>k -i, x jtk -i)), 
we introduce yi ik = Xi,k ~ x i,k-i — &&{.fi{pi,k-i) + Pi,j9i x i,k-i, x j,k-i)}i an d expand using the Taylor expansion. 



h(x it k) exp | ^fi( x i,k) + ^2p'i,j9i x i,k,Xj,k)\ 



{xi,k+i - Xi,k) + O(At) + 



E h ™y™k- 



(28) 



Therefore, we obtain 



/Qp r 

dxi. k h(x hk )- | e =o = C j dyi t k^2h m y^ k e^p(-y^ k /2a 2 At) 

°£i,k J 

= c J dy^Y^i^Aty^K^eM-y 



i,k) 



(29) 
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Since the value of this integral becomes if m is odd, the leading order of this integral is obtained from the term 
?7i = 2, and has a mag nitude of 0{At 3 / 2 ). Since the magnitude of this contribution is smaller than O(At), we can 
neglect this term at the limit At — > 0. 

Therefore, if ki, k2, and k$ are unequal, then the value of the integral of the second term in Eq. (|20|) becomes 
0. In the case where k\ = hx — k%, d 3 P(e)/de 3 is negligible as At — ► from a similar argument. In general, the 
contribution from the term d m P(e)/de m is negligible for odd m, because it is expressed in the form e~ Sl (xi t k — 



j Xj,k—i 



))) x (nonsingular function). 



Similarly, we can prove that [d m /dei 1 ,k 1 ■ ■ ' d e i m ,k m ]P is nonzero if and only if each e^k appears 2n times in the 
delimiter. From these results, we conclude that the contribution of the term I = 3 to Z is of the order pf ^ . In addition, 
we conclude that the contribution from the term I = 2m + 1 is of the order p\ ■. The contribution from the term 
I = 2 m (2m' + 1) is also estimated by the Stratonovich-Hubbard transformation. In the case of the random network 
model, the correction due to these terms is small. For example, the contribution from the term / = 6 to the MSR 
generating functional term is estimated using 



exp 




iAtXi t kg(xi,k-i, x^k- 



I 180 



drij exp 




-1> Xj,k-l) 



(30) 

The dispersion of r^j is y/piJ/QO and the contribution of this term is much smaller compared to that from the 
Gaussian fluctuation obtained from the term 1 = 2. 

From the discussion based on the Stratonovich-Hubbard transformation and perturbation, we demonstrated that 
the MSR generating functional for the dynamics of a random sparse network model is almost identical to that for the 
dynamics of a random Gaussian network. In the next section, we demonstrate that the above analysis is consistent 
with the result of a numerical simulation of the Kuramoto transition in a network model. 



IV. EXAMPLE: THE KURAMOTO TRANSITION 



In the previous section, we developed a general scheme to approximate the dynamics of the random sparse network 
and found its dynamics can be described by 

N 

x\ = fi(xi) + y^(pi,j + rij)g(xi,Xj) + (31) 

3 

when pi j « 1. In this case, the distribution of r^j is provided by a Gaussian with dispersion a = ^/Pij. In this 
section, we apply this approximation to the dynamics of oscillators in random networks. 
We consider a random network of oscillators 

ddi/dt = Lu t + K^2 a %o a[n ( ] - e i) ( 32 ) 

3 

where 9i and u>i represent the phase and velocity of the oscillator i, respectively. The value of dj j is 1 if nodes i and j 
are connected, otherwise it is 0. We note that the random Gaussian matrix needs to be considered as symmetric. We 
consider the case where the distribution of u>i is given by g(ui) = (1 / \/2iw\J) exp((w — ujq) 2 / 20^) , and a^j represents 
a random network with its mean degree ko. The above discussion suggests that the dynamics of this network can be 
approximated using the following equation: 

ddi/dt = u>i + ^ (j± + Ksmie,- - 9i), (33) 

j 

where the distribution of Tij is given by P(r^j) = ^/N2irko exp(—Nrfj/ (2fco)). This model is similar to the dynamic 
glass model proposed by Daido|2l|. However, the mean interaction between oscillators is positive in our model, while 
it is in Daido's model. 

It is difficult to calculate analytically the dynamics of this globally coupled model. In this section, we present 
the numerical results for the random sparse and random Gaussian networks. For this simulation, we set N = 1000, 
2<r5 = 1.0, ojq = 0, and k = 10. The result obtained is averaged over 50 different networks. 

First we examine the coupling dependence of the order parameter r. In our previous paper, we defined the order 
parameter r as r = 'Yl li kie l9i / ~^2,ki for a random sparse network model. However, it is difficult to define such an 
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FIG. 1: The coupling dependence of the order parameter for random sparse and random Gaussian networks. 

order parameter for a random Gaussian model. In this paper, we therefore use (e lS ) as the order parameter for a 
random Gaussian network, and r = Yli^i e / 12^i f° r a random sparse network model. Although these two order 
parameters are distinct, the difference between them is small, because the distribution of degree has a strong peak 
at k = ko for a random sparse network. The values of r are plotted in Fig. ^ f° r both the networks for the range 
K = 0.02-0.20. In both these models, the order parameter remains almost constant for K values less than 0.1. There 
is a rapid increase in r for K values greater than 0.1. The values of r coincide qualitatively for these two models. 
When K = K c , a sharp transition, given by r oc \]K — K c , is observed in the mean- field approximation. This sharp 
transition gets smeared out in a random sparse network. The Gaussian model approximates this smearing well. 

The order parameters being identical is not unusual, because their obtained values were close to those obtained 
using the mean- field theory. We now explain the distribution of velocity d6i/dt. In the mean-field approximation, 
d6/dt has a J-function-like peak at d6/dt = 0. However, if the coupling between the oscillators is random, the strong 
peak at dd/dt will get smeared. In Fig. [3 the distribution of d6/dt for sparse random and random Gaussian networks 
is plotted. At K = 0.02, there is an absence of synchronization and the distribution of d6/dt is Gaussian-like. On 
the other hand, the oscillators are well synchronized and the distribution has a strong peak at dO/dt = when 
K = 0.16. For the present study, we focus on the distribution at K = 0.10. This value of K is close to the critical 
point, and we suggest that the large fluctuation appears at this point. In the case of sparse networks, the peak at 
dO/dt = is sharper at k = 0.10 than at k = 0.02. The same tendency is observed in the case of a Gaussian network. 
For example, we observe that P{— 0.1 < d9/dt < 0.1) = 0.137 for a sparse random network. This value is close to 
P(— 0.1 < d6/dt < 0.1) = 0.131 obtained from a random Gaussian network. This consistency in the observed value 
suggests that a random sparse network can be approximated by a Gaussian random network. 

Finally, we present the distribution of the phase 9 for both networks. The phase distribution in the (cj, 9) plane at 
K = 0.16 is shown in Fig. |3| Although the coupling strength is sufficiently large for synchronization, the phase does 
not entirely lie on a single line obtaind from the mean-field approximation method, 9 = arcsin(cj/.fi* k$r). In order 
to observe the dispersion around the mean-field line, we present the phase distribution of oscillators with < 0.05 
in FigQ] In this region, | axcsin(ui/Kkor)\ is less than 0.05 and the u> dependence of the phase distribution can be 
neglected. In both the models, the phase distribution lies in a wide range of 9. The dispersion a for these two figures 
is a 2 = 0.105 and 0.122 for the random sparse and the random Gaussian network, respectively. Since these two values 
coincide qualitatively, random Gaussian network is a good approximation of the random sparse network. 

V. CONCLUSION AND DISCUSSION 

In this paper, we studied the dynamics of a random network model using the path-integral approach. We identified 
that the mean-field approximation is the lowest-order approximation of p^j and I = 1, as shown in Eq. Ijllfl . We also 
demonstrated that the contribution of the term I — 2™ can be described by the fluctuation of coupling in the globally 
coupled approximation method. The contribution of the odd I terms is difficult to estimate, though it is of the order 
pf j. We applied these general results to the Kuramoto transition, and observed a good agreement with numerical 
simulations. 

The path-integral approach developed through this study is a general one and is applicable to dynamics of any 
random network. In particular, if the precise result for a randomly coupled model is known, a good approximation 
can be obtained for random sparse network models. There are several models, such as the replicator model|23|. for 
which the exact results are known for a Gaussian random network . Our analysis proves that the dynamics of random 
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FIG. 2: The distribution of d6/dt at K = 0.02,0.10, and 0.16 for a random sparse network (upper) and a random Gaussian 
network (lower). 





FIG. 3: Phase distribution in sparse and Gaussian random networks at K = 0.16. 



sparse networks can be easily obtained for such models. 

The analysis presented in this study is limited to the dynamics in a random network model. In the case of another 
network model, we need to include the higher-order terms to evaluate the MSR generating functional. It is usually 
difficult to carry out such a calculation. However, our result provide much informations regarding the validity of 
the mean-field approximation. For example, the mean-field approximation is applicable if PijPk,i <C Pij <C 1. On 
the other hand, such an approximation is not applicable to the dynamics of a highly clustered network. In such a 
network, PijPj,kPk,i ~ 0(pi^pj.k), PijPk,i <C Pi.j cannot be assumed and the contribution from the neglected terms 
needs to be calculated. It is usually believed that the dynamics of networks with high clustering coefficients cannot 
be approximated using the mean-field approximation method because of the high clustering coefficient. However, 
our analysis reveals that the validity of the mean-field approximation methods depends on the value of PijPk,i and 
Pi j. For example, the mean-field approximation method cannot be applied to the square-lattice model even if the 
clustering coefficient is zero, because the value of PijPk,i can be as large as that of pij. 

We also discuss other studies conducted on the Kuramoto transition in random network models. Restrepo et al. 
examined the mean-field theory and studied the Kuramoto transition jlfjj . They concluded that synchronization 
occurs when K satisfies the relation K > 2/irg(Q)\, where A is the largest eigenvalue of the network matrix a^j and 
<?(0) is the density of the oscillators at u> = 0. They stated that the mean-field approximation, which was developed by 
us in previous papers, functions only when cx fej, where r$ is the local field defined as rj = (y]j <hj^ ^^^ti where 
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FIG. 4: Phase distribution in random sparse and random Gaussian networks at K — 0.16. 



(• • -) t means the average over a long time interval. However, they did not explain the reason why this assumption is 
valid some random network models, though they stated that there exists some relationship between the eigenvectors 
of ciij and degree of each node. In this paper, we demonstrated that the mean-field theory is an approximation that 
considers only the term I = 1 in the MSR generating functional. In this case, the mean-field approximation coincides 
with the discussion obtained from the largest eigenvalues, because the largest eigenvalue of the matrix pi^ — XiXj is 
Y). x 2 and its eigenvector v is given by v — {x\,X2, • • -x n ). In the random network model, pij = kikj/N(k), where 
ki and kj are the degrees of the nodes i and j. Therefore, the largest eigenvalue of this matrix is (k 2 )/(k), and 
the critical condition for synchronization in the mean-field approximation becomes identical to that in the discussion 
based on eigenvalues. In order to examine the applicability of the mean-field approximation, the term I = 2 should be 
considered. In the case of a random matrix, the largest eigenvalue with a dispersion p is expressed by 2yJNp — 2 -J (k) 
based on Wigner's semicircle lawj^. This result suggests that the mean-field approximation can be applied if 
yj (k) <C (k 2 )/(k). In order to examine this, we consider the matrix M + G, where M is the matrix obtained from 
the mean-field approximation and G the Gaussian random matrix, i.e., the distribution of each element of the matrix 
is Gaussian with dispersion y/k. As observed earlier, the largest eigenvector v of matrix M satisfies the condition 
Mv = \v, where A = (fc 2 ) / (k) . On the other hand, \Gv\ is of the order \/2fc|?;|, because all eigenvalues of G lie between 
— v2fc and y/2k. Therefore, \(M + G)v\ equals approximately X\v\ based on the assumption that A ^> y/k, and the 
direction of (M + G)v is approximately identical to w. It should be noted that all vectors u that are perpendicular 
to v, i.e., (u,v) = 0, satisfy the condition Mu = 0. This implies that \(M + G)u\ < y/2k\u\ -C A|m|. Therefore, the 
largest eigenvalue and corresponding eigenvector of M + G can be approximated as (k 2 ) / (k) and v — (k\, ■ ■ ■ ,k n ), 
respectively. Therefore, the mean-field approximation is a suitable approximation if (k) is sufficiently large. In the case 
of a scale-free network, the spectrum density differs from that suggested by Wigner's law and the above-mentioned 
conclusion should be modified. However, this discussion suggests that the validity of the mean-field approximation is 
determined by the largest eigenvalues of the mean-field matrix M and random matrix G. If the largest eigenvalue of 
matrix G is as large as A, the mean-field approximation is not valid. Based on the idea presented in this paper, the 
claim made by Restrepo et al. implies that the term I = 2 must be included in order to discuss the critical behavior 
more accurately, especially in the case of a scale- free network with 7 = 2. Therefore, their work was not a denial of 
the mean-field theory, but an extension of it. 
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